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自 组织 多 目标 粒子 群 优化 算法 
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摘 要 : 针对 多 目标 粒子 群 优化 算法 收敛 性 和 多 样 性 难以 平衡 的 问题 ， 提 出 一 种 利用 问题 的 结构 信息 来 解决 多 目标 问 
题 的 自 组 织 多 目标 粒子 群 算法 。 通 过 自 组 织 映 射 网 络 发 现 种 群 和 非 支配 解 集 分 布 的 结构 ,构造 出 当前 粒子 的 邻 域 关 系 ， 
从 邻 域 中 选 出 非 支配 解 ， 从 而 引导 种 群 局 部 和 全 局 的 搜索 。 提 出 了 精英 学 习 策略 ， 通 过 对 精英 粒子 进行 变异 ， 引 导 算 
法 跳出 局 部 最 优 。 实 验 结 果 表 明 ， 所 提 算 法 可 以 兼顾 收敛 性 和 多 样 性 ， 有 效 地 解决 多 目标 优化 问题 。 
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Abstract: In order to keep the balance between convergence and diversity of multi-objective particle Swarm optimization 
algorithm, the paper proposed a self-organizing multi-objective particle Swarm algorithm (SMPSO), which utilized structure of 
population to solve multi-objective optimization problems. The distribution of the population and non-dominated solutions 
founded by self-organizing map network helped to construct the particle neighborhood relations , so as to improve the 
algorithm’s local and global search ability by selecting non-dominated solution as leader in the neighborhood. The elite learning 
strategy could help population jump out local optimum by doing mutation on elite particles. The experimental results show its 
ability to keep convergence and diversity, and show its effectiveness to solve the multi-objective optimization problems. 
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0 ”引言 性 算法 ， 一 类 是 进化 算法 。 进 化 算法 是 一 类 以 种 群 为 基础 的 带 
有 随机 搜索 方式 的 方法 ， 具 有 经 过 不 断 演化 而 获得 将 近 真 实 帕 

在 现实 应 用 中 ， 有 很 多 问题 需要 同时 优化 多 个 目标 ， 两 个 ” 累 托 前 沿 的 能 力 ， 因 此 被 广泛 应 用 到 多 目标 问题 的 求解 上 ， 如 
及 以 上 目标 的 问题 叫做 多 目标 优化 问题 (multi-objective ”NSGA-II “中 六 到 引用 浙 ， 多 目标 差分 (DE) “看 引 用 渐 、 多 目标 粒子 
optimization problems,MOP)。 这 些 目标 之 间 往 往 是 相互 冲突 的 ， ”和 群 (MOPSO) 算 法 9 在 以 种 群 为 基础 的 算法 中 , 粒子 群 
一 个 目标 的 性 能 提高 会 导致 其 他 目标 性 能 的 降低 ， 因 此 同时 优 算法 具有 快速 收敛 的 特点 二 于,。 对 于 多 目标 粒子 群 优化 
化 多 个 目标 是 很 困难 的 。 不 同 于 音 问题 ， 多 目标 问题 。 算法 ,平衡 收 敛 性 和 多 样 性 是 至 关 重 要 的 。 因 此 提出 了 很 多 改 
得 到 的 不 是 一 个 解 ， 而 是 由 许多 非 支 配 解 组 成 的 解 集 ， 称 为 由 ” 进 策略 :随机 选择 策略 ， 从 外 部 存储 中 随机 选择 一 个 全 局 最 优 
累 托 解 集 。 帕 累 托 解 集 在 目标 空间 对 应 的 目标 值 便 组 成 了 帕 累 ”gbest 是 简单 且 高 效 的 策略 ; 超 体积 贡献 ,拥挤 距离 等 也 可 以 作 
托 前 沿 。 由 于 多 目标 问题 得 到 的 解 可 能 有 无 限 多 个 ， 所 以 找 出 ”为 选择 gbest 的 指标 ; 一 个 种 群 分 不 同 偏好 的 子 种 群生 
所 有 的 解 是 不 实际 的 。 在 实际 操作 中 ， 找 出 一 定数 量 的 非 支配 。” “， 一 个 负责 全 局 搜索 ， 男 一 个 负责 局 部 收敛 ; 文献 错误 ! 未 找 
解 , 尽 可 能 履 盖 并 靠近 真正 的 帕 累 托 前 沿 , 且 这 些 解 分 布 均匀 ， 到 引用 源 。 中 粒子 不 同 维度 向 不 同 的 个 体 学 习 ， 可 以 增加 种 群 
是 本 文 所 需要 解决 的 挑战 性 难题 。 多 样 性 。EPSO “ee 引入 了 均衡 因子 ， 以 提高 种 群 全 局 搜 
解决 多 目标 问题 的 方法 有 两 大 类 ， 一 类 是 传统 的 数学 解析 索 能 力 。AGE-MOPSO “fen 利 用 种 群 收敛 状态 ， 自 适应 
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调整 惯性 权重 和 学 习 因 子 的 值 , 以 达到 探索 和 开发 的 最 佳 平衡 。 
但 这 些 策 略 很 少 利 用 到 多 目标 优化 问题 的 特性 ， 即 在 一 定 的 条 
件 下 ，m 个 目标 的 多 目标 问题 在 目标 空间 的 帕 累 托 前 沿 及 决策 
空间 的 帕 累 托 解 集会 形成 m-1 维 的 片段 连续 的 流 型 形状 外 。 基 
于 多 目标 问题 的 这 种 规则 性 , 文献 错误 ! 未 找到 引用 源 。 提 出 将 


法 few。 它 通过 模仿 岛 类 运动 规则 在 搜索 空间 中 寻找 最 
多 位 置 。 与 其 他 进化 算法 相似 ， 搜 索 由 种 群 个 体 来 组 成 ， 这 些 
搜索 个 体 叫做 粒子 。 粒 子 由 两 个 特征 即位 置 和 速度 信息 组 成 。 
假设 当前 种 群 由 图 个 粒子 组 成 ,粒子 的 搜索 范围 是 维 超 空间 。 
粒子 ie(,2,.N) 的 位 置 和 速度 信息 分 别 表示 为 


吾 
Xi (Ki Warr Wn) 和 Vv = (Wa Va Vy) 区 入 更 新 公式 如 下 : 


自 组 织 映射 网 络 应 用 于 种 群 的 重组 操作 中 ， 有 效 地 提高 了 多 


标 算法 的 性 能 。 但 目前 尚未 有 人 将 自 组 织 映射 网 络 与 多 目标 粒 
子 群 优化 算法 结合 。 本 文 提 出 的 自 组 织 引 导 机 制 的 自 组 织 多 目 
标 粒 子 群 优化 算法 (self-organizing multi-objective particle swarm 
optimization algorithm，SMPSO)， 利 用 自 组 织 映射 网 络 可 以 将 
高 维 信息 映射 到 低 维 网 络 ， 且 保持 特征 分 布 不 变性 的 特点 ， 映 
射 当 前 种 群 和 帕 累 托 解 集 的 结构 信息 ， 为 粒子 构建 新 的 邻 域 天 
系 ， 引 导 粒 子 的 搜索 过 程 ， 以 此 平衡 多 
性 和 多 样 性 。 


1 ”相关 背景 知识 


1.1 多 目标 优化 问题 


标 粒 子 群 算法 的 收敛 


一 个 连续 的 多 目标 优化 问题 可 以 表述 为 如 下 形式 : 
Min F(X) =(f(),f.00,.. ,00) 人 


St. = ) EO 
其 中 : Q=[4,,b 了 是 搜索 空间 的 可 行 域 ,x 是 决策 空间 9 中 的 
n 维 决策 变量 ; m 是 待 优 化 目标 的 个 数 。 函 数 F:Q 一 R" 定义 
了 m 个 待 优化 目标 函数 由 决策 空间 Q 到 目标 空间 R" 的 映射 。 


些 多 目标 问题 的 相关 定义 如 下 : 
定义 1 帕 累 托 支 配 性 。 决 策 变量 x* 支配 决策 变量 y( 记 为 


x<y ) 当 且 仅 当 满足 以 下 条 件 : 
(Vie{l 2 m): £0 < £0) 
人 (3j et{l, 2 mm}: fi(X) < £7)) 
即 在 m 个 目标 函数 中 ,x 的 每 个 目标 函数 值 都 不 大 于 y 的 目标 
函数 值 ， 且 x 至 少 有 一 个 目标 函数 值 小 于 y 的 目标 函数 值 。 
定义 2 帕 累 托 解 。 决 策 变 量 x 是 帕 累 托 解 : 
yeQ:y<x (3) 
即 在 可 行 域 空间 内 ， 不 存在 任何 一 个 粒子 支配 粒子 x。 
定义 3 帕 累 托 解 集 (Pareto optimal set，PS) 定义 为 
PS={xeQ| -HyeQ:y< (4) 
定义 4 帕 累 托 前 沿 (Pareto ffont，PF) 定义 为 
PF ={F()|xePS} (5) 
1.2 多 目标 问题 的 规则 性 
在 一 定 条 件 下 ， 一 个 连续 的 m 个 
形成 m-1 维 的 流 型 ， 如 0 所 示 。 
鉴于 这 种 规则 性 ， 有 很 多 构造 分 布 模型 的 算法 二 下 
例如 用 主 成 分 分 析 法 构造 模型 和 下 9"， 利 用 模型 产生 新 解 ， 
但 这 种 构造 模型 的 方法 ， 效 果 依 赖 于 模型 参数 的 选取 。 
1.3 ”粒子 群 优化 算法 
粒子 群 优化 算法 是 基于 种 群 的 随机 搜索 的 全 局 寻 优 进 化 算 


(2) 
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标 问题 的 PF 和 PS 会 


V(t+1)=w*yv,(t) + cn (x 


pbess, 一 五 CD)+C2D Xow 一石 (OD) (6) 
X(t+1)=x (CD)+V LT+D (7) 


其 中 : t 是 迭代 代数 ; w 是 可 以 用 来 平衡 局 部 和 全 局 搜索 的 惯性 
因子 , 较 高 的 惯性 因子 给 先前 经 验 大 的 权 值 而 有 益 于 全 局 搜索 ， 
反之 , 偏向 于 局 部 搜索 ，c, 2 代表 两 个 加 速 因子 影响 种 群 的 收 
敛 速 度 ; ,是 在 (0,1) 正 态 分 布 之 间 产 生 的 随机 量 ; pbest 和 
gbest 分 别 是 个 体 最 优 和 全 局 最 优 解 。 
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图 1 多 目标 问题 的 规则 性 


多 目标 粒子 群 算法 中 ,寻找 的 不 是 唯一 的 gpest， 而 是 用 邻 
域 最 优 的 个 体 nbest 代替 gpest。 式 〈6) 改 为 如 下 形式 : 


V(t+D) =w*y,(t)+ Ci (Xppest, x (D) + Npes, —xX(1)) (8) 


2 ” 自 组 织 多 目标 粒子 群 优化 算法 SMPSO 


本 文 提出 的 SMPSO 算法 ， 主 要 是 利用 自 组 织 映射 网 络 
(self-organizing mapping network, SOM ) 发 现 种 群 和 多 目标 问 
题 帕 累 托 前 沿 的 结构 信息 ， 然 后 引导 粒子 的 飞行 。 同 时 结合 
精英 学 习 策略 和 多 项 式 变异 , 增加 种 群 多 样 性 。 SMPSO 流程 如 
0 所 示 。 


(初始 化 ) 
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图 2 SMPSO 的 流程 


2.1 SOM 模型 为 PSO 粒子 构建 邻 域 信息 
神经 网 络 与 进化 算法 的 结合 大 多 是 用 进化 算法 优化 神经 网 
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络 的 参数 守 吕 下 效用 本 文中 利用 自 组 中 映射 网 络 为 粒子 群 算 
法 中 的 粒子 构建 邻 域 关 系 ， 从 而 引导 粒子 的 飞行 。 

如 0 所 示 , SOM 包括 一 个 输入 层 和 一 个 表示 层 。 假设 输入 
数据 是 n 维 ， 表示 层 有 入 = Ni x, 个 神经 元 。 每 个 神经 元 
ue (1,2,…, 入) 赋予 位 置 z* = (zy, 芝 ) 和 权 值 w= (wi ,wh,…, Ww) 。 


有 表示 层 
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图 3 SOM 的 拓扑 结构 示意 图 


SOM 的 流程 分 为 两 个 阶段 : a) 学 习 阶 段 ， 随 机 选择 训练 数 
据 ， 根 据 欧 氏 距离 选择 获胜 神经 元 ， 更 新 获胜 神经 元 及 其 邻 域 
神经 元 的 权 值 ，b) 聚 类 阶段 ， 将 测试 数据 映射 到 神经 元 ， 相 似 
的 数据 会 映射 到 相 邻 神经 元 。 

构建 的 SOM 模型 ， 神 经 元 的 权 值 w = (Cs) 的 初 
始 值 与 种 群 粒 子 的 位 置信 息 相 同 ， 每 迭代 一 代 更 新 一 次 。 进 化 
过 程 中 得 到 的 新 的 非 支 配 解 的 位 置信 息 是 SOM 的 训练 数据 
(TS), SOM 更 新 之 后 , 将 种 群 (POP) 和 外 部 存储 集 (EXA) 


的 粒子 位 置 映射 在 网 络 结构 上 ， 根 据 神 经 元 的 位 置信 息 ， 构 建 
粒子 的 邻 域 ， 从 而 选择 合适 的 引导 粒子 abest。 有 具体 步骤 如 算法 
1 所 示 。 

2.2 算法 1 nbest 的 选择 


输入 : POP，EXA，TS， 
习 半 径 ou 。 
输出 : ww = (5) ， 每 个 粒子 的 nbest_i。 


权 值 w= 的 0) ,学 习 率 11 ,学 


a) 计算 每 个 神经 元 之 间 的 距离 


| 三 和 


，, 根据 距离 排序 ， 


经 元 u 的 邻 域 神经 元 集合 了 ， 是 距离 神经 元 u 第 k" 近 的 
神经 元 。 


“ 


梁 静 ， 等 织 多 目标 粒子 群 优化 算法 
: 1 一 CO *(] 。 = 其 三 | 
其 中 : 7(D) = *( TREE ; GOOD)=aos( 1) ; MAXgen 是 最 


大 迭代 代数 。 

d) 将 待 分 类 数据 POP、EXA 输入 ， 根 据 式 “9) 判断 获胜 神 
经 元 。ind2unit 记录 POP 每 个 粒子 的 获胜 神经 元 , unit2nbest 
记录 每 个 神经 元 上 对 应 的 nbest。 

e)nbest i 记录 粒子 x 的 nbest 备 选 库 

for 1=1:NP 

u= ind2unit (i); /粒子 i 的 获胜 神经 元 u 

nbest_i= unit2nbest {了 };/ 神 经 元 u 的 邻 域 神经 元 上 对 
应 的 nbest / 

end 

如 0 所 示 ， 当 前 粒子 i 的 获胜 神经 元 是 w， 根 据 神经 元 的 
位 置 确 定 邻 域 神经 元 ， 邻 域 神经 元 上 对 应 的 EXA 中 的 个 体 则 
组 成 nbest 的 备 选 库 nbest_i。 参 照 快速 排序 法 ， 根 据 占 优 性 和 
拥挤 距离 ， 选 择 排 序 靠 前 且 稀 疏 的 个 体 为 粒子 i 的 nbest。 相 较 
于 将 目标 空间 划分 成 几 部 分 的 分 解 策略 竺 下 并 局 汪 。 氏 未 坑 到 引用 浙 。 
此 策略 有 以 下 优点 : 

a)SMPSO 选 出 来 的 粒子 在 决策 空间 距离 粒子 较 近 , 产生 
效 子 代 的 效率 更 高 。 

b)SMPSO 在 决策 空间 为 粒子 构建 邻 域 , 整个 种 群 有 不 同 的 
引导 粒子 ， 保 持 了 种 群 的 多 样 性 。 

cjSOM 的 结构 是 根据 目前 得 到 的 非 支配 解 集 不 断 学 习 的 ， 
权 值 信息 不 断 更 改 ， 更 靠近 PF， 产 生 的 解 更 符合 问题 的 结构 。 
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妈 4 ”npbest 的 选择 


2.3 ”扰动 变异 策略 

对 于 多 目标 粒子 群 算法 ， 保 留 非 支配 解 是 至 关 重 要 的 。 因 
为 它们 是 当前 搜索 到 的 最 佳 位 置 ， 具 有 一 定 的 引导 意义 。 但 | 
于 缺乏 对 待 优化 问题 的 先 验 知识 , 当 种 群 处 于 收敛 停滞 状态 四 


Py 


b) 选择 训练 数据 X= (4,…%,…%,) ETS ， 判 断 获 胜 神经 元 : 


CT 
w=arg min|—w|, (9) 


其 中 ， 上 | 是 欧 氏 距离 。 


c) 更 新 获胜 神经 元 及 其 邻 域 神经 元 wsU) 的 权 值 : 


U={u|ll<u<N ^|z -2 <oln)} 


(10) 


wa =w +nD*expCle’ zw) 


启 


有 可 能 是 陷入 了 局 部 最 优 。 本 文中 除了 使 用 外 部 存储 集 之 外 还 
引入 了 精英 学 习 策略 ， 在 目前 得 到 的 最 优 位 置 处 进行 变异 ， 有 
利于 全 局 最 优 解 的 进一步 寻找 。 

在 为 粒子 x 选择 nbest 之 后 ， 进 行 子 代 更 新 ， 更 新 之 后 的 
粒子 按照 一 定 概率 进行 精英 学 习 。 
j=rand(1,n)<pr; 


x(j) =nbest(j) + normmd(0, pr”) x (b,-a)) (11) 


其 中 : 0 和 D 是 粒子 x 第 闫 维 的 上 下 界 ; pr 的 值 随 着 迭代 


次 数 的 增加 ， 从 0. 2 线性 减少 为 0. 05。 在 进化 初期 ， 较 大 的 
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变异 范围 有 利于 算法 的 全 局 搜索 ， 随 着 和 从 代 的 进行 ， 算 法 需要 的 快速 排序 法 对 备 选 库 的 粒子 进行 排序 ， 选 择 排序 靠 前 且 稀 下 


在 局 部 范围 找到 更 精确 的 解 ， 因 此 小 的 变异 范围 更 适合 。 正 态 。 的 粒子 作为 nbest。 
分 布 的 概率 密度 如 图 5 所 示 。 由 0 可 知 ， 当 pr 值 较 大 时 ， 产生 (0) 子 代 更 新 。 
的 随机 值 较 大 的 概率 比较 大 ， 从 而 粒子 的 变异 范围 较 大 ， 有 利 在 SMPSO 中 ， 粒 子 群 更 新 公式 为 
于 产生 较 大 的 扰动 ， 随 着 迭代 的 进行 ，pr 值 减 小， 产生 的 随机 v=wxv+randxc x(pbest ~ +randxc,x(nbest—*) (14) 
数 在 均值 0 附近 的 概率 会 越 来 越 大 ， 变 异 范围 减 小 ， 此 时 的 ee. 
nbest 已 接 近 于 真实 的 最 优 位 置 ， 小 范围 的 变异 更 有 利于 寻找 YY MXsn 
精确 解 。 速度 边界 处 理 : 超过 边界 的 值 赋予 边界 值 。 
8- 无 三 X 二 TY (15) 
ER 之 后 ， 按 照 一 定 概率 执行 精英 学 习 策略 ， 参 考 式 (11)。 
6 po1 | 位 置 边界 处 理 ， 超 过 边界 的 值 赋予 边界 值 且 速度 置 反 。 
‘pr=0.05 进行 多 项 式 变异 操作 ， 参 考 式 (12) 和 (13)。 
4 A | (d) 更 新 EXA 和 TS。 


外 部 存储 集 EXA 的 更 新 策略 ， 根 据 占 优 性 和 拥挤 距离 排 
序 ， 保 留 排序 靠 前 且 稀 玻 的 粒子 ; 
训练 数据 TS=unique(EXA,TS), 即 只 把 新 产生 的 加 入 EXA 

的 粒子 作为 训练 数据 。 
图 5 正 态 分 布 的 概率 密度 氏 dd) 结束， 否则 返回 步 又 o)。 
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同时 ， 为 了 避免 种 群 陷 入 局 部 最 优 ， 引 入 了 多 项 式 变 异 。 En 
if randO<l/n 3.1 实验 设置 
x(D)=x(D+6)x(b)—a)) (12) 3.1.1 测试 函数 
a tg dd 本 文采 用 文献 错误 ! 未 找到 引用 源 。 中 的 GLT1-GLT6 和 文 
= 1 3) 献 错误 ! 未 找到 引用 源 。 中 DLTZ2、DLTZ4 和 DLTZ5 为 测试 函 
0 0 数 。0 列 出 了 测试 函数 的 决策 变量 和 Pareto 前 沿 的 信息 。 
其 中 : rnd 是 在 〈0,1) 间 的 随机 数 ，eta 是 变异 分 布 指数 。 表 1 用 于 算法 性 能 比较 的 多 目标 优化 测试 问题 
2.4 SMPSO 算法 的 具体 步骤 测试 问题 ”目标 维 数 ”变量 维 数 变量 范围 Pareto 前 沿 特点 
在 SMPSO 中 ， 利 用 SOM 模型 为 粒子 构建 邻 域 ，SOM 权 GLTI1 2 10 [0.H]x[-11] 汪 线性 不 连续 
值 随 着 种 群 每 迭代 一 次 进行 一 次 更 新 ， 训 练 数据 只 采用 新 加 入 GLT2 2 10 [0,]]x[-1,1J”* ”上 是 PF 连续 
外 部 存储 集 EXA 的 粒子 ， 节 约 计算 量 且 避免 数据 训练 过 度 。 GLT3 2 10 [0,1] x[-1,1]"™ 凸 PF 连续 
在 子 代 更 新 中 加 入 了 扰动 策略 ， 平 衡 了 收敛 与 多 样 性 。 对 粒子 GLT4 2 10 [0.1]x[-11] 和 一 凸 PF 不 连续 
位 置 进行 了 边界 处 理 ， 边 界 粒子 速度 置 反 可 以 避免 粒子 集中 于 GLTS 3 10 [0,1 了 x[-b41T 了 ”点 PF 连续 
边界 位 置 。 算 法 2 给 出 了 SMPSO 算法 的 具体 步骤 。 GLT6 10 [0,1 x[-1,1J”” ”点 PF 不 连续 
2.5 算法 2 SMPSO 的 具体 步骤 DTLZ2 3 12 [0,17 叫 形 ， 连 续 
输入 : 种 群 大 小 NP，PSO 参数 w cl, c2. DTLZ4 3 位 [0,17 凹 形 非 均 匀 
输出 : 外 部 存储 集 EXA。 DTLZ5 3 12 [0,]] 线 型 连续 


a) 初 始 化 种 群 和 SOM 网 络 权 值 。 
初始 化 种 群 NP 个 粒子 POP={x,x,.…,x*}, EXA=POP; SOM 3.1.2 比较 算法 


I SOD OP 为 了 比较 算法 (SMPSO) 的 性 能 ,选取 三 种 代表 性 的 多 目标 
0% 优化 的 对 比 算法 。 这 些 算 法 中 有 NSGAIT 和 *! 吉 到 引用 要。、 
b) 根 据 神经 元 之 间 的 距离 此 一 z"|| ,计算 神经 元 u 的 邻 域 ee 
错误 ! 未 找到 引用 源 。 和 MOEA/IGD-NS 错误 ! 未 找到 引用 源 。 各 算法 的 实验 参 > 数 都 
神经 元 集合 了 ， 了 是 距离 神经 元 u 第 ke 近 的 神经 元 。 按照 对 应 参考 文献 设置 。 
oO) 循环 。 3.1.3 实验 参数 
(a)SOM 权 值 更 新 : 根据 TS 及 式 (9)(10) 更 新 网 络 权 值 。 实验 中 GLT1-GLT4 是 2 目标 的 多 目标 问题 ， 种 群 大 小 设 


(b) 参 考 算法 1， 为 粒子 选择 nbest 备 选 库 ， 按 照 NSGA-II 为 100, SOM 的 网 络 结构 设置 为 一 维 (1x100 ), 外 部 存储 容量 


均 为 100， 测 试 函数 的 最 大 评价 次 数 为 50 000。GLT5-GLT6 和 1 DLZ4-DLZ5 的 参考 点 7 = (2,2,2) 。 

DTLZ 系列 是 3 目标 的 多 目标 问题 ， 种 群 大 小 设 为 150，SOM 实验 中 JGD 和 “HV 的 值 都 用 来 衡量 种 群 的 收敛 性 和 分 布 
的 网 络 结构 设置 为 二 维 (15x10 )， 外 部 存储 容量 均 为 150, 测 。 性 ，1GD 值 越 小 (HV 的 值 越 大 )， 则 得 到 的 PF 越 好 。 得 到 的 
试 函数 的 最 大 评价 次 数 为 70 000。 每 种 算法 在 所 有 测试 函数 ” 1GD 值 很 小 且 HV 值 很 大 ， 则 得 到 PF 非常 接近 真实 的 PF， 且 


站 


上 均 独 立 运行 25 次 。 能 完全 覆盖 真正 的 PF。 
3.1.4 性 能 指标 3.2 ”实验 结果 及 分 析 
为 了 评估 算法 获得 的 近似 PF 的 收敛 性 和 均匀 性 ， 本 文采 SMPSO 算法 中 需要 更 新 拓扑 结构 的 权 值 , 且 在 一 定 概率 下 


用 反 转 世代 距离 (inverted generational distance，IGD) 和 超 体积 ”进行 变异 。 为 了 比较 各 种 算法 的 时 间 复 杂 度 ， 将 6 个 比较 算法 
值 (hyper-volume，HV) 作 为 性 能 评估 指标 。 令 了 为 真实 Pareto ”在 9 个 多 目标 测试 问题 上 独立 运行 25 次 的 耗 时 均值 列 于 0 内 。 


最 优 解 集 , P* 是 得 到 的 Pareto 最 优 解 集 , 则 7GD 的 计算 如 下 : 可 以 看 出 SMPSO 算法 的 时 间 复 杂 度 比 NSGAII 和 IM-MOEA 
0O 》 村 司 一 数量 级 ， 示 的 测 i 问题 本 | 
JGD(P',P) = 2* seP'd0,P) a 的 大 ， 但 保持 在 同一 数量 级 ， 在 3 目标 的 测试 问题 上 耗 时 比 

四 2LBMOPSO，pccsAMOPSO 和 MOEA/IGD-NS 相对 较 少 。 


0 和 0 列 出 了 6 种 比较 算法 在 9 个 多 目标 测试 问题 上 的 
IGD 和 HV 性 能 指标 值 。 其 中 , 平均 值 (mean) 和 标准 差 值 (stq) 是 


其 中 ，d(x',P) 是 x 到 P 中 任何 一 个 点 的 最 小 距离 ，|P | 是 P* 


的 基数 。 同一 算法 在 同一 测试 问题 上 独立 运行 25 次 的 统计 结果 ; 同时， 
HV 的 计算 方式 为 采用 粗 体 字 和 灰色 背景 表示 所 有 对 比 算法 在 每 一 个 测试 问题 中 


的 最 小 1GD 值 和 最 大 HV 值 。 每 种 算法 对 同一 个 问题 的 排名 及 
综合 排名 列 在 表 内 。 同 时 用 秩 和 检验 ,5% 的 显著 性 水 平 比较 各 
村 其中; = Ci 六 ) 是 目标 空间 里 被 任何 帕 累 托 前 沿 占 优 的 参考 。 算法 与 SMPSO 的 显著 性 差异 。 分 别 用 “1”、“§ ”和 “” 表 
点 ，VOL() 是 勒 贝 格 测度 。GLT1、GLT3 的 参考 点 r=(2,2) ， 。 示 算法 SMPSO 比 对 比 算法 较 好 , 较 差 和 没有 明显 差异 的 情况 。 
< GLT2 的 r= (2.1D) ,GLT4 的 参考 点 7 = (2,3) , GLT5-GLT6、DLZ2 


HV(P,r)=VOL JLAGO, NIX-[f Dn) a7) 


表 2 各 种 算法 在 每 个 测试 问题 上 的 平均 耗 时 /s 


2LBMOPSO NSGAII pecsAMOPSO IM-MOEA MOEA/IGD-NS SMPSO 
GLT1 62.14 12.38 82.72 18.49 5.49 46.3 
GLT2 63.75 10.86 98.07 19.57 5.94 45.35 
SC GLT3 68.88 人 97 143.67 24.22 24.38 45.36 
GLT4 72.42 12.03 84.86 20.41 10.78 50.41 
一 GLTS 138.1 19.91 138.20 25.68 183.45 76.11 
GLT6 140.84 20.99 602.35 23.66 176.82 74.56 
DTLZ2 134.24 20.16 635.679 20.12 426.54 76.41 
DTLZ4 131,33 20.7 543.43 20.45 37261 79.87 
DTLZ5 130.42 20.38 605.28 22.92 390.28 78.56 


表 3 测试 函数 性 能 度量 1GD 的 均值 (标准 差 ) 


函数 2LBMOPSO NSGAII pccsAMOPSO IM-MOEA MOEA/GD-NS SMPSO 
mean 1.17E-01'[6] 4.35E-02 [3] 8.12E-02'[4] 3.61E-02 [2] 9.88E-02'[5 2.38E-02[1] 
GLTI 
std 4.40E-02 3.07E-02 2.41E-02 1.95E-02 5.22E-03 2.25E-02 
mean 1.12E-01'[3] 4.28E-02'[2 7.12E-01'[5] 6.59E-01'[4] 1.95E+00'[6 3.82E-02[1] 
GLT2 
std 2.20E-01 2.89E-03 4.67E-01 2.22E-01 1.60E-01 1.88E-03 
mean 6.64E-02'[5] 2.81E-02'[2 1.12E-01'[6] 4.17E-02'[3] 6.53E-02'[4 9.32E-03[1] 
GLT3 
std 2.17E-02 1.46E-02 4.46E-02 8.31E-03 5.53E-03 5.32E-03 
mean 1.32E-01'[4] 1.16E-01'[3 2.27E-01'[5] 6.14E-02'[2] 2.62E-01'[6 4.17E-02[1] 
GLT4 
std 7.04E-02 7.80E-02 9.81E-02 4.16E-02 1.50E-02 4.90E-02 
mean 5.14E-02'[4] 6.51E-027[5 8.88E-02'[6] 4.86E-02 [3] 3.43E-02°[1] 4.77E-02[2] 
GLTS 


std 2.70E-03 4.57E-03 1.84E-02 2.96E-03 1.25E-03 2.57E-03 
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mean 4.99E-02'[4] 5.83E-02'[5] 7.34E-02'[6] C3.18E-02™[2] 2.71E-02 [1] 4.12E-02[3] 
GLT6 
std 1.77E-03 5.37E-03 3.59E-02 1.26E-03 1.56E-03 3.34E-03 
mean 6.15E-02'[4] 1.00E-01'[6] 5.27E-02° [2] 7.78E-02'[5] 4.29E-02°[1] 5.51E-02[3] 
DTLZ2 
std 2.46E-03 7.48E-03 1.69E-03 3.50E-03 3.65E-04 2.04E-03 
mean 6.28E-02'[2] 3.15E-01' [6] 1.05E-01'[4] 6.73E-02' [3] 2.12E-01'[5] 5.40E-02[1] 
DTLZ4 
std 1.64E-03 7.99E-02 1.16E-01 2.38E-03 3.02E-01 2.05E-03 
mean 4.36E-03'[5] 3.69E-03 [4] 3.37E-03°[2] 1.57E-02'[6] 2.80E-03°[1] 3.68E-03[3] 
DTLZ5 
std 2.68E-04 5.14E-04 2.16E-04 1.96E-03 2.19E-04 2.41E-04 
Rank 4.111 4.000 4.444 3.333 3.333 1.778 
1 § /> 9/0/0 7/0/2 7/0/2 6/0/3 5/3/1 
表 4 测试 函数 性 能 度量 HV 的 均值 (标准 差 ) 
函数 2LBMOPSO NSGAII pccsAMOPSO IM-MOEA MOEA/IGD-NS SMPSO 
mean 一 一 [6] 3.26E+00 [3] 3.17E+00'[4 3.28E+007[2] 3.17E+00'[5] 3.31E+00[1] 
GLTI1 
std 一 一 6.84E-02 5.38E-02 3.44E-02 2.42E-02 5.30E-02 
mean 1.84E+01'[4] 1.97E+01 [2 1.83E+017[5 1.90E+01[3 1.64E+01'[6] 1.98E+01[1] 
GLT2 
std 1.11E-02 7.03E-03 6.47E-0 2.34E-01 5.51E-01 1.01E-02 
mean 3.93E+00'[5] 3.94E+00'[4 3.91E+00'[6 3.94E+00'[3 3.94E+00 [2] 3.95E+00[1] 
GLT3 
std 4.03E-03 2.75E-03 3.70E-02 1.74E-03 1.13E-03 1.06E-03 
| 二 
A 入 mean 4.86E+00'[2] 4.85E+00'[3 4.40E+00'[6 4.81E+00'[4: 4.42E+00'[5] 4.97E+00[1] 
GLT4 
std 1.27E-01 1.83E-0 9.14E-0 1.35E-01 2.07E-02 2.59E-02 
mean 3.88E+00'[3] 3.87E+00'[4 3.85E+00'[5 3.81E+00'[6 3.88E+00'[2] 3.90E+00[1] 
GLTS 
std 1.31E-02 2.35E-02 2.75E-02 2.03E-02 1.56E-02 1.56E-02 
mean 3.85E+00'[4] 3.84E+00'[6 3.86E+00" [2 3.85E+00'[5 3.89E+00 [1] 3.86E+00[3] 
GLT6 
std 1.08E-02 1.42E-02 3.32E-02 2.34E-02 1.27E-02 1.37E-02 
mean 3.48E+00'[4] 3.47E+00'[5 3.49E+00 [2 3.48E+007[3 3.46E+00'[6] 3.49E+00[1] 
DTLZ2 
std 2.46E-03 7.48E-03 1.69E-03 4.25E-02 2.94E-02 3.59E-03 
mean 3.48E+00'[5] 3.57E+00°[1] 3.55E+00’[2 3.52E+00 [3] 3.47E+00'[6] 3.50E+00[4] 
) DTLZ4 
€ std 1.64E-03 7.99E-02 1.16E-01 2.63E-02 2.08E-01 2.05E-03 
mean 2.84E+00 [1] 2.84E+00'[5] 2.84E+00" [2 2.84E+00'[4] 2.83E+00'[6] 2.84E+00[3] 
DTLZ5 
std 2.68E-04 5.14E-04 2.16E-04 2.30E-02 6.19E-03 2.16E-04 
Rank 3.778 3.667 3.778 3.667 4.333 1.778 
TS 8/0/1 7/1/1 5/1/3 7/1/1 7/1/1 


从 0 中 可 以 看 出 ，SMPSO 算法 对 大 多 测试 问题 都 有 很 好 ”试问 题 上 效果 都 比较 好 , 如 在 DTLZ4 上 效果 较 差 。 综 合 排 名 和 


的 性 能 度量 。 对 于 测试 问题 GLT1-GLT4 和 DTLZ4，SMPSO 入 显著 性 水 平 显示 ，SMPSO 对 大 部 分 测试 问题 都 有 较 好 的 多 样 
法 得 到 的 1GD 均值 最 小 。 对 于 测试 问题 GLT5, SMPSO 的 IGD 性 。 

值 仅 大 于 MOEA/IGD-NS 算法 ,对 于 GLT6、DLTZ2 和 DLTZ5， 为 了 进一步 显示 各 种 算法 求 得 最 优 解 集 的 收敛 和 多 样 性 的 
SMPSO 算法 的 IGD 指标 排名 第 三 。 表 中 的 Rank 值 是 算法 对 所 情况 ,以 测试 问题 GLT1、GLT3、GLT6、DTLZ4-DTLZ5 为 例 ， 


有 测试 问题 的 平均 排名 ， 可 以 看 出 SMPSO 的 排名 最 高 。 且 秩 ”将 测试 函数 运行 25 次 结果 ， 按 1GD 值 排名 中 间 的 一 次 运行 结 
和 检验 的 结果 表明 ，SMPSO 在 大 多 测试 问题 上 明显 优 于 其 他 ， 果 显 示 如 0 所 示 。 


对 比 算法 。 GLT1 GLT1 
0 中 “~~” 表 示 2LBMOPSO 算法 对 测试 问题 GLT1 的 结 avopsol [ee | 

果 较 差 ， 性 能 指标 HV 出 现 负 值 ， 不 作 统计 。SMPSO 算法 对 于 和 “| 

GLT1-GLT5 和 DTLZ2，HV 值 均 为 最 大 。 对 于 GLT6，DLTZ5 ~ | 全 本 

测试 问题 ，HV 均值 排名 第 三 。 当 然 ，SMPSO 不 可 能 在 所 有 测 ls 
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(c) GLT6 的 PF 对比 
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6 测试 函数 PF 对 比 
从 0 中 可 以 看 出 , 无 论 是 2 目标 还 是 3 目标 问题 , SMPSO 
得 到 的 PF 均 可 以 较 好 地 覆盖 真正 的 PF， 相 对 于 对 比 算法 ， 


SMPSO 能 较 好 的 处 理 


| 


息 ， 且 网 络 结构 的 整 


E PF 不 连续 的 情况 ， 如 GLT1 和 GLT6， 


在 SMPSO 中 的 SOM 网 络 结构 可 以 较 好 地 发 现 问 题 的 结构 信 


本性 有 利于 种 群 的 全 局 搜索 ， 因 此 找到 的 


PF 没有 缺失 。GLT3 虽然 PF 是 连续 的 , 但 是 测试 函数 的 复杂 ， 


导致 找到 的 PF 能 完 


DTLZ4 的 问题 上 , SMPSO 得 到 的 PF 虽然 能 较 全 面 


全 履 盖 真正 的 帕 累 托 前 沿 有 一 定 的 难度 。 


地 履 盖 真正 


的 PF， 但 图 中 显示 分 布 的 均匀 性 略 差 。 DTLZ5 测试 问题 比较 


简单 ， 大 部 分 算法 得 
为 了 验证 


到 的 帕 累 托 前 沿 能 够 靠近 真正 的 PF。 


SOM 构造 邻 域 关 系 的 合理 性 , 0 将 以 GLT2 问 


ChinaXiv 会 作 期 刊 | 

未 用 杂 前， 等 由 央 吕 多 则 村 了 可 许 全 办 法 
题 为 例 , 列 出 在 进化 过 程 中 神经 元 上 对 应 的 粒子 及 粒子 的 nbest a 
备 选 库 在 决策 空间 和 目标 空间 的 分 布 。 al | 一 一 

从 图 中 可 以 看 出 , 在 进化 的 过 程 中 , 为 粒子 选择 的 nbest 不 lk A 用 
论 在 PF 或 是 PS 空间 ， 都 是 外 部 存储 集 EXA 中 距离 粒子 个 体 at 人 em 
比较 近 的 ， 有 利于 局 部 搜索 ， 增 加 产生 高 质量 粒子 的 可 能 性 。 | 
粒子 a 和 b 对 应 的 神经 元 在 拓扑 结构 中 距离 较 近 ， 它 们 映射 的 0 
粒子 距离 也 比较 近 。 同 理 ，c 粒子 对 应 的 神经 元 在 拓扑 结构 上 ， | : | 
距离 较 远 ， 则 a、b 距离 c 粒子 及 其 邻 域 较 远 。 这 也 是 SOM 的 ie 
特性 之 一 ， 即 相似 的 个 体 聚 在 一 起 。 每 进化 一 代 之 后 ， 权 值 学 = 
习 一 次 , 因此 权 值 的 调整 也 是 慢 慢 累积 的 。 图 (b) 中 粒子 c 在 PF Pr 
空间 分 布 距 离 邻 域 较 远 , 是 因为 SOM 的 映射 关系 是 在 PS 空间 图 8 不 同 pr 初始 值 对 各 个 测试 问题 IGD 值 的 影响 


建立 的 ， 所 以 反映 在 PF 空间 会 有 所 偏差 。 


SMPSO 中 根据 多 目标 优化 问题 的 规则 性 ，2 目标 问题 , 采 
和 pF Bg 用 一 维 (1x100 )SOM, 3 目标 问题 , 采用 二 维 (15x10 )SOM。 
二 人 丰 了 、 
1 i 外 为 了 验证 SOM 拓扑 结构 对 算法 性 能 指标 IGD 和 HV 的 影响 ， 
8 1 % 将 SOM 的 拓扑 结构 作出 如 下 改变 : 2 目标 问题 改 用 二 维 
a ol 1 和 (10x10 ) SOM，3 目标 问题 改 用 一 维 ( 1x150 ) SOM， 将 算 
-0.5. 本 a 、 
4 < i \ 法 独立 运行 25 次 的 平均 值 结果 显示 如 0 所 示 。 
、 lh . a 
i r | 
i \ Gi | 1 0.08 E 
0 一 人 2 0 上 一 维 SOM 
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fl 
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2 aeons 2 Sp (a) SOM 拓扑 结构 对 1GD 的 影响 


[可 一 维 SOM| 
(第 400 代 时 粒子 的 分 布 | 
， 是 当前 外 部 存储 集 EXA; 人 
妇 是 粒子 a， 第 40 号 神经 元 对 应 的 粒子 ， 去 是 粒子 a 的 npest 备 选 ， a 
国足 粒子 b， 第 45 号 神经 元 对 应 的 粒子 ， 〇 是 粒子 b 的 nbest 备 选 ; EE 
国 是 粒子 <， 第 80 号 神经 元 对 应 的 粒子 ，[_] 是 粒子 c 的 nbest 备 选 。 | | | 
图 7 不同 粒 子 的 邻 域 分 布 中 | | | | | | | ] 
SMPSO 算法 中 精英 变异 策略 中 的 参数 pr 初始 值 的 选取 ， 9 GL GLT2 GLT3 GLT4 GLTS GLTE DTL22 DTLZ4 DTLZS 


测试 问题 


在 参考 了 数学 模型 的 概率 分 布 之 外 ， 另 外 进行 了 实验 的 验证 。 
如 0 所 示 , 设 定 pr 最 小 值 为 0.05, pr 初始 值 分 别 设 为 0.1、02、 
0.3、0.4 和 0.5， 折 线 图 显示 没有 某 一 个 pr 初始 值 可 以 使 所 有 
测试 问题 的 1GD 值 均 为 最 小 , 但 是 pr 初始 值 设 为 0.2 时 可 以 使 
大 多 测试 问题 的 1GD 值 较 小 , 且 变 异 概率 较 小 可 以 降低 计算 复 
杂 度 。 因 此 pr 初始 值 设 为 0.2, pr 线性 减 小 到 最 小 值 为 0.05 而 
不 是 0， 是 为 了 在 算法 的 后 期 仍 能 产生 一 定 的 扰动 。 


(b) SOM 拓扑 结构 对 HV 的 影响 


图 9 SOM 拓扑 结构 对 IGD 和 BV 的 影响 
从 图 9 (a) 中 可 以 看 出 ， 对 于 2 目标 的 测试 问题 GLT1、 

GLT3，1GD 值 在 二 维 SOM 的 情况 下 较 小 ，GLT2、GLT4, 在 一 
维 SOM 的 情况 下 较 小 ，3 目标 的 测试 问题 在 二 维 SOM 时 IJGD 
值 较 小 ， 从 (b〉 中 可 以 看 出 ，SOM 拓扑 结构 的 改变 对 2 目标 
问题 的 HV 值 影响 不 大 ，3 目标 问题 的 HV 值 在 二 维 SOM 的 情 
况 下 略 高 。 因 此 ， 总 的 来 说 ，2 目标 问题 对 一 维 或 二 维 拓扑 结 
构 的 SOM 无 明显 偏好 ,3 目标 问题 较 偏 好 二 维 拓 扑 结构 的 SOM 。 
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4 ”结束 语 


通过 结合 多 目标 问题 的 规则 性 和 自 组 织 映 射 网 络 的 映射 关 
系 不 变性 ， 本 文 提出 了 一 种 新 算法 。 该 算法 为 多 目标 粒子 群 算 
法 构建 了 一 种 新 的 邻 域 关系 ， 为 每 个 粒子 选取 了 更 合适 的 引导 
粒子 ,同时 加 入 了 精英 学 习 机 制 ,算法 兼顾 了 收敛 性 和 多 样 性 ， 
对 大 多 的 多 目标 问题 都 能 得 到 分 布 均匀 且 靠 近 真正 PF 的 解 集 。 
在 后 续 研 究 中 ， 将 该 算法 应 用 到 高 维 多 目 标 问 题 上 ， 并 进一步 


改进 算法 的 性 能 。 


tot 
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